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Abstract 

The effect of the most general Z — > 37 vertex in the reaction e + e~ — > 3j is studied with a 
particular attention to LEP searches. We give exact analytical expressions including realistic 
cuts for the signal and present a detailed analysis based on a Monte Carlo that includes the 
effect of the irreducible 37 QED cross section. As special applications we discuss the effect of 
heavy scalars, fermions and gauge bosons and comment on the "monopole" connection. 
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1 Introduction 



The second phase of LEP will soon be in operation. The achievements of the first phase 
have already been formidable as we have learnt a great deal from the precision measure- 
ments at the Z peak on the parameters of the Standard Model, SM . Some of the future 
analyses will need to exploit the full statistics to be accumulated at the end of the LEP1 
runs| in order to put stringent constraints on rare processes that are usually the hallmark 
of New Physics. One such process is the decay of the Z into three photons. Within 
the SM this process is so rare that even the once-discussed High- Luminosity option of 
LEP1 f| would not be able to see it. Indeed, branching ratio for this process calculated 
within the SM is Br sm (Z — > 37) ~ 5.4 1CT 10 0. This makes it an ideal candidate to 
search for contributions that are beyond the SM . It is interesting to note that this kind 
of quartic vertex, Z37, is an example of an "anomalous" boson self-coupling which can 
be probed at LEP1 whereas the second phase of LEP will be mostly dedicated to the 
tri-linear WWj, WWZ couplings. The study of this coupling thus provides a nice bridge, 
when moving from LEP1 to LEP2, in the general topic of the direc^ investigations on the 
vector boson self-couplings. From the experimental point of view, another motivation for 
studying this decay is the observation that the description of this rare process is, as we 
shall see, identical to the description of the scattering of light by light parameterised by a 
4-7 coupling. Although QED has been with us for a long time, this non-linear effect has 
not been directly experimentally investigated. As far as we are aware, one has only exper- 
imentally studied Delbriick scattering and photon splitting near a heavy nucleus at low 
energy [[|. One should therefore not miss any opportunity of investigating a very similar 
kind of physics, especially since the theoretical situation at LEP1 with the Z — > 37 is 
much cleaner. Thus one should take full advantage of the resonance enhancement at the 
Z peak which provides an alternative to study New Physics contributions to an Abelian 
non- linear effect. 

Before the advent of LEP, it had been argued |5], || that this decay could check the 
hypothesis of a composite Z. The argument heavily borrows from the p system and vector 
meson dominance 0. To explain the universality of the weak coupling in this picture, one 
adapts the vector dominance idea to the W triplet by implementing a W°— 7 transition ||. 

^By the end of 1994 each LEP collaboration has collected about 140 pb _1 (~ 5.5 10 6 Z) of data 
samples but only a fraction of these has been used in published analyses of rare processes. At the end of 
the LEP1 runs one hopes to accumulate about 200 pb _1 per experiment Q. 

*Of course, LEP1 is of an unsurpassable precision when it comes to the investigation of the vector 
bosons self-energies and hence on what might be considered as bi-linear anomalous couplings. 
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The universality of the QED electromagnetic coupling is thus transmitted to the W/Z- 
fermion coupling. The difference between the p system and the W is that the strength 
of this transition, directly related to the weak mixing angle sw, is not small as would 
be expected for an electromagnetic transition but of order 1/2. Since this is large, one 
should expect that other electromagnetic transitions be enhanced 0. More specifically, 
in this picture, one can build up on the "strength" of the Zj transition to predict a large 
decay of the Z into three photons. Of course, the Z can not decay into two on-shell 
photons (Yang's theorem f[0| ) and in any case the Z2*y vertex violates C and V and thus 
could be forbidden. However, Z — > 37 is perfectly allowed. Arguing along these lines, a 
non-relativistic bound state calculation || of the Z — > 37 leads to a branching ratio of 
the order of 10~ 5 . This is in fact the present order of magnitude the LEP experiments 



HI 12, 13, 14 1 have set, as a limit, on this decay. 



Recently De Puijula |1| has suggested that this decay may be large if it is induced 
by "monopoles". This idea has already been put forth [[1(J a few years ago to motivate 
a non- negligible 4-7 coupling. The aim of this paper is, first, to give the most general 
framework for the study of the Z —>■ 3j where the above two examples ( so-called "com- 
posite" Z and "monopoles" ) will be seen to be specific cases of the most general effective 
Lagrangian. We will then give a detailed analysis on the signature of the Z — > 37 decays. 
For the most general Lagrangian we give exact analytical formulae for various distribu- 
tions including the effect of the correlation with the initial e + e~state. We also derive 
analytical formulae taking into account experimental cuts that can be used for a quick 
estimate of the acceptances of the 4 LEP detectors. Finally we conduct a detailed dis- 
cussion based on a Monte Carlo analysis which includes not only the signal but also the 
irreducible e + e~ —>■ 37 purely QED background as well as the effect of the interference. 
We show which distributions are most sensitive to the presence of the Z — * 37 coupling. 
Finally a short discussion about the effect of this coupling far away from the Z resonance 
is presented. In the Appendix we collect some simple and compact expressions for the 
standard and non-standard helicity amplitudes based on the technique of the inner spinor 
product that were implemented in our matrix element Monte Carlo generator. 



2 Effective Lagrangians and Models for Z — > 37 

Probably the oldest example of an effective Lagrangian is the celebrated Euler Lagrangian 
fT7| that describes the self-interaction of photons. The first non-trivial interaction de- 
scribes a 4-7 vertex in its leading part and accounts for the scattering of light-by-light. 
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The basic idea of the effective Lagrangian is to give a general description of a phenomenon 
even if one does not know its origin or the underlying theory behind it. All we need to 
construct the ensuing operators is the known symmetries at low energies which are, in 
this example, contained in the U(l) gauge invariance. For the 4-7 vertex one is led to 

C t}f. = ^71 (HF^n 2 + HF,uFn 2 ) with = V Q/3 F Q/3 . (2.1) 



We remark that even though one is dealing with an effective Lagrangian, and restricting 
oneself to the leading operators, the Lagrangian is very constrained since it contains 
only two operators. These are necessarily of dimension eight. The mass M is directly 
related to the scale of New Physics. Non-leading operators contain extra derivatives that 
give correction factors of order (E^/M) 2 to the dominant terms we are considering, and 
therefore their effect is very negligible. We can also argue that if these terms were to be 
taken into account then this would also mean that multi-photon amplitudes 67,.... etc 
would not be negligible and thus the New Physics would be more conspicuous. 

Different models, or rather types of heavy particles, give definite predictions for (3 2 ,z- 
For instance, considering the case of a massive particle (mass M) of unit charge that 



contributes to the one-loop 4-7 amplitude, we have in the case of fermions |L7| [[! 



P( = i/3f = -1. (2.2) 



For a charged scalar the coefficients /3j's are different with |L8 



$ = 7/3| = — — . (2.3) 

H2 H 3 16go V J 

For completeness, we can also give the contribution of a heavy charged gauge vector boson. 
Here as the spin increases the values of the coefficients increase as well. This is not only 
due to the counting of degrees of freedom. Although the coefficients for the f3's in the 
case of scalars and fermions have been known for a long time Hl8| , it is only in the last 



few years that the issue of the vector, more specifically the W^-loop, has been investigated 
19| . The f3's can be inferred from [ IS ] : 



gv = qv = (24) 

H2 160 3 160' v ' 1 

It is instructive to realise that if instead of considering ^2,3 one takes the combination 

f3± = (3 2 ±(3 3 (2.5) 
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the counting of the degrees of freedom, with the correct Bose-Fermi statistics factor, is 
reflected only in /3_: 



3 5 1 -c 1 



P+ 180 ' P ~ 240 
#=^j ; P F . = -2(3 S _ (2.6) 

We will have more to say about the extraction of these coefficients and about the so-called 



-2p 3 )(F„ u F^)(F„ v Zn + 4p 3 (Z^F vX F Xp FW)) . (2.7) 



super symmetric relations |2(J between the /3_ elsewhere pT|j . 

The Z — > 37 vertex is constructed along the same mould. In fact, since we are replacing 
only one photon by a Z the general leading-order effective Lagrangian has exactly the 
same structure. Therefore, we choose to parameterise the Z37 anomalous coupling in the 
form 

4o\ 
M 4 V 

The factor of 4 in the definition of the Z37 coupling compared to the case of the 47 
coupling is such that the two Lagrangians lead to the same strength, and expression in 
momentum space, for the Z37 and 47 coupling modulo the overall change /3j — ► 

We refrain here from giving specific examples of the coefficient ^2,3 for the different spin 
species as, beside the electric charge of the particles, we need to know their hypercharge 
or SU{2) weak assignment. Nonetheless the relations between (3 2 and /3 3 are maintained 
in the same form as in the all 4-7 case (Eq. |2.2[ - 12.3| ), apart from the case of the spin-1, 
where there is a subtlety. For the latter one needs to know both the quantum number 
of the pure gauge part and the Goldstone part of the massive vector bosons. In practice 
the pure gauge part is by far dominant. Since these are rather technical observations and 
calculations we leave this discussion to a forthcoming paper |2l] . 



Again it is worth working, instead of the two independent constants $ 2 an d $3, with 
the combination 

P± = fa±h- (2-8) 

The reason is that /3 + and /3_ contribute to different helicity amplitudes and therefore 
they do not interfere. It is easy to see the reason why this is so. For simplicity it is 
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best to stick with the 4-7 Lagrangian and rewrite it in terms of the self-dual, F^ u , and 
anti-self-dual, F~, fields which for photons correspond to definite helicity states, namely 

F± u = ~(F^ ± i F" v ) ; F_- F + = F» v F + ^ = 0. (2.9) 

In this case one has 

C%, = ^ (2/3 + FlFl + /?_ ((Fir + (Fir)) • (2.10) 

It is then obvious that the terms in (3 + contribute to different helicity amplitudes than 
those with /?_: formally taking a configuration with only either a self-dual or an anti-self- 
dual (either F + —>■ or F~ — > 0) the (3 + does not contribute. For the same reason, in the 
process Z — > 37 based on the identical Lagrangian, only /?_ contributes to the amplitude 
where all 37 have the same helicity. 

For the discussion we will consider 5 typical one-parameter "models" of (3: 

• (32,3'- As stressed above, since it is only the combination (3± that appears at the 
helicity amplitude level, it is not possible to disentangle a model with only /3 2 from 
a model with only (3 3 in three-photon production. The model with (3 2 has, in an 
abuse of terms, been referred to as a "composite" Z. j3 2 could, more appropriately, 
be associated with the effect of a heavy scalar having a two-7 and Z'j coupling, 



while P3 with that of a pseudo-scalar. Two of us [22| have devoted a preliminary 



analysis to the effect of the operator fi 2 some time ago. 
a model with (3 + but /3_ = 
a model with /?_ but f3 + = 

a model that we will refer to as the (3^ model with (3 + = —11/3 This is the same 
model as that considered in HTsR. It corresponds to the effect of a heavy fermion 



which in [|J| is considered to be a monopole. In [JT5J] the exotic fermion is assumed 
to be a SU(2) singlet. With this information it is straightforward to relate the 47 
and Z37 strengths as 

ks = — Pis (4, = 1-40 (2.11) 

c w 



where /3<f 3 are given by Eq. |272|. Taking the fermion to be a point-like monopole 
(free isolated magnetic charge) one is also led to assume that the strength of its 
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coupling to the photon, g m , is strong. In fact Dirac's quantification condition 
relates this coupling to the electric coupling through 

eg m = 27m (n = integer). (2-12) 



Therefore in our parameterisation one has to allow for the change e — > g m in Eq. \2.1\ . 
Although one may be tempted by giving a monopole motivation to this type of 
coupling, it remains that one still has two parameters: the mass of the monopole 
and the magnetic strength, g m (or quantification number n) which are lumped in 
one effective parameter M/n, exactly as what one has with the effective Lagrangian 
(Eq. |2.7| ). The discriminating relation Eq. [2.11| indicates the spin- 1/2 nature of the 



point-like monopole. Note, for the record, that the idea of a monopole to induce 
a large 4-7 coupling is not new; it has been exploited some time ago by Ginzburg 
and Panfil [jl(J . They also considered the case of a spin-0 monopole which again can 
be accommodated by the general effective Lagrangian through the discriminating 
relation = ^-^23 ( an d e — ► g m ); however, from what we remarked earlier, for 
the same quantification condition the spin- 1/2 monopole gives a larger contribution. 

In fact, by far, the most interesting scenario is that of the effect of a heavy spin-1 
boson. This could be a model with j3 + = 28/3_. This corresponds to the effect of 
a very heavy vector boson whose pure gauge component and Goldstone component 
have the same quantum numbers. This may be considered as the vector boson 
version of the "monopole" considered by De Rujula ||15|| . In practice this model 



is very similar to a model with (3- = in view of the much larger value of /3+, 
especially if one keeps in mind that the main contributions of these operators, as we 
shall see, are quadratic in the couplings. For a thorough discussion on these models 
and others we refer to |2T . 



In the monopole vein, one should not dismiss the eventuality that due to the neces- 
sarily very strong coupling of the monopoles, bound states may form with a mass that is 
much smaller than their combined constituent masses. One interesting possibility would 
be that the ground state is of spin-0 that decays predominantly into 77 and Zj, there- 
fore dynamically realising models with ^2,3 . In the remaining of the paper we will steer 
away from considerations pertaining to the manifestation of the monopoles. This is not 
only because one would have to argue that their masses are, against what is generally 
but not universally believed, much below the GUT scale (in order to evade a Lilliputian 
unobservable effective Z — > 37 coupling), but also because one has to quantify their 
effect on the two-point functions that are extremely well measured at LEP1 and contrast 
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the information with the one given by Z — > 37 |. We will stick to the general unbiased 
approach of studying the manifestation of a general Z — > 37 coupling referring to the 
models as defined through the relation between j3 + and (3-. We take the view that the 
"monopole" is just a possible paradigm for the existence of a not too small Z — > 37 
coupling. 



We also note for further discussion that the authors of reference |25| have also given a 
parameterisation in terms of two operators but have only studied the double integrated 
Z — > 37 partial width and the totally integrated width. To make contact with their 
parameters we note that their constants G\p are related to ours as 

Gl = ^(A-A) = -^ G^ h . (, 13) 
3 Analytical Formulae for Y[Z — > 37) and e + e~ — > 37 

From the general Lagrangian, it is an easy task to compute the Z width into three photons 
as well as the correlation with the initial state through the cross section e + e~ — > 37. We 
will start by providing analytical formulae for the partial width and will give compact 
analytical formulae including realistic experimental cuts, that could be considered as 
canonical cuts for this process. The latter can be used to quickly estimate the acceptances. 



Figure 1: The t-channel QED process and the s-channel anomalous Z — > 37 contribution 
to three-photon production. 




From the start, we would like to point out, before exhibiting any formula, that since 

the effective operator is of dimension 8, constructed out of field strengths, the s-channel 

§ Moreover, beside the controversial "lightness" (~ TeV) and point-like nature of the Dirac monopole, 
consistency of QED in the presence of a Dirac monopole has been questioned (23 . 
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produced photons, see Fig. 1, will be hard photons that will not be collinear to the beam. 
This is in contrast to the t-channel QED background (see Fig. 1). This is a general 
characteristics of photons that probe New Physics (See P^| ). In the case at hand, we 
should therefore expect that the distribution in the least energetic photon, for instance, 
will be more sensitive to the presence of the Z — > 37 vertex. 

Specialising to the decay Z — > 37 and taking as variables the reduced photon energies 
Xi = 2Ei/M z , the double distribution in the energies writes: 

1 dr 



r dxidx 2 
120 



2 + /3 2 _)(l-x 1 )(l-x 2 )(l-x 3 ) + 



3/? 2 + 5/3 2 

01 + 2© ((1 - x0 2 (l - x 2 f + (1 - xx) 2 (l - x 3 ) 2 + (1 - x 2 ) 2 (l - x 3 ) 2 ) } 
60 w {01 + 2/3!) ((1 - Xl fx\ + (1 - Xl fx\ + ((1 - * 3 ) 2 * 2 ) 



3(3 2 + + 5/3 2 

- 2 / 3 2 (l-x 1 )(l-x 2 )(l-x 3 )}, (3.1) 
where the totally integrated partial width is calculated to be 



r^^L^YfWim. (3. 2) 

18 7T 3 \ M 2 J \ 120 / K ' 

We disagree with the expression given in [15] in the case of the fermion monopole. 



although we confirm the expression for the normalised differential cross section. We find 
a rate that is 12 times smaller than the one given in |l5| |"] . On the other hand, we 



confirm the direct one-loop calculation of the supersymmetric charged Higgses given in 
|26| . Needless to say that for this weakly interacting heavy particle the effect is hopelessly 
beyond observability. We also agree with the expression given in |25 in terms of Gi,2- 



We note, in connection with the monopole case, that had we considered a spin-1 rather 
than spin- 1/2 exotic, the spin-1 gives a factor of about 120 enhancement (if one keeps the 
same other quantum numbers)! 

From these expressions we see that present experimental limits [|T^, [Tj|, Br(Z — > 37) ~ 
10~ 5 place a bound on the strength of the Z37 of the order of 

4(^Y/W4. (3.3) 



^ This would mean that the extracted L3 [|T3) limit on the mass of the monopole should be corrected 
by a factor ~ 2. Moreover from the MC that we will provide it will be possible to implement the correct 
acceptance factor for this model. 
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Therefore, one can easily argue that a conventional weakly interacting particle (ft ~ 1) 
can not be probed through the Z — > 37 decay. This is not to say that we should not 
look for these decays; this would be like arguing that lacking any model that gives large 
WW^j/Z anomalous coupling one should not probe the WW^/Z couplings at LEP2! 

The single-photon energy distribution is 



ldr_ 20 ^( ( ~o2 , 52 w-, ^ , P+ + 2 P 
f dx ~ 3ftl + 5ft 



■x 



(ftl + ft 2 _)(l -x) + ^- (20 - 40a; + 21a; 2 ) . (3.4) 



It is more telling and appropriate to distinguish the photons by ordering them according 
to their energy; we will label them as E\ > E 2 > E 3 (Ei is the energy). The energy 
distribution of the least energetic photon is 

1 dr 12 
Tdx~ 3 3 ftl + 5ft 2 _ 

{xl [(ft 2 + + 2ft 2 _)(l + 3x3(1 - x 3 ) + 24(1 - x 3 ) 2 ) - 501(1 - x 3 )] 9(^ - x 3 ) 

+ (2 - 3a; 3 ) [(ft 2 + + 2ft 2 _)(l + 3(1 - x 3 )(l + (1 - x 3 )(A - 23x 3 + 27a; 2 )) 

5ft 2 _(l - x 3 )(l - 3(1 - x 3 ) 2 )] 9(2/3 - x 3 )6(x 3 - 1/2)} , (3.5) 

where 6(x) is the step function. 

The average energy of the least energetic photon, E 3 , is 

M z 1 /751~ 2 367 ~ 2 



This shows that independently of ft-,ft+ the average energy of the photon is almost 
the same for the two independent operators, i.e., §.22\M Z for ft + and 0.218 for ft_. As 
expected these are hard photons as compared to the expected average energy of the softest 
photon one gets for a typical QED 37 event. Therefore, it is relatively easy to disentangle 
the QED and Z — > 37 events; however it is not easy to distinguish between different 
types of New Physics. This will be made clearer when we show the results of the full 
simulations including cuts. It is also difficult to differentiate between the models on the 
basis of the average energy of the most energetic photon and the "medium" one. 

It is relatively straightforward to implement the cuts within the event plane as done 

in the LEP experiments. Requiring the separation angle between any two photons (6ij) 

to be above a certain value and requiring a cut on the energy of the photons, such that 

. 2E- 

9ij > 2VArcsin 5 ; Xi = > e, (3.7) 

we calculate the acceptance to be 
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^— - 1 ^ ~ 1 - f + 2/3 2 ) — + (/3 2 + / 3 2 ) ^ + (4/3 2 + 5/3 2 ) — 

r + 5/52 7 VP+ 7 v H+ 11 7 

+ {5(31 + 9/3!) 15£ 4 



(3.8) 

We refrain from giving the exact expression, not only because it is lengthy and unin- 
spiring but also because the typical values of the cuts applied by the LEP experiments 
are such that the above approximate expression is more than enough. Note for instance 
that the cut on the energy is of order e 4 reflecting once again the preferentially energetic 
photon characteristics of the New Physics and the quartic dependence of the differential 
width on the energy. This is the reason we have kept terms up to fourth order in the 
cut-offs. At this order there is "no-mixing" between the cuts on the opening angles and 
the energy. With E { > M z /20 and % > 20° (See OPAL 0), one does not lose more 
than a few percent of the signal. The message then is that one can afford applying harder 
cuts than the above without much affecting the signal; even with cuts that are twice as 
large as the typical LEP experiments cuts our approximate formula, Eq. \i.8[ is excellent. 
The main purpose of the "mild" LEP1 experimental cuts is to bring about a decent and 
clean 37 sample (which in the absence of any new physics is purely QED). With higher 
statistics we could afford cutting harder. It is worth emphasising that the acceptance is 
sensibly the same for different manifestations of New Physics, i.e., whether one takes /3 + 
or /3_ . Again it has to do with the fact that the dependence on the photon energies for 
different choices of (3 are not terribly different. 

Introducing the correlations with the initial state, we have been able to obtain some 
exact results including cuts. For the sake of generality we also give the cross section for 
polarised electron beams. 

For the 37 final state we have that 

z 9l 

v( e+e L,R 37) = 2 2 L ' R 2 (Tunp. (3.9) 

9l + 9r 

where <T unp . is the unpolarised cross section and 

9l = ~^ + Sw ; 9r = s 2 w . (3.10) 

We introduce, for a centre-of-mass energy, yfs, the Z propagator via 

(s- M 2 7 ) - iT z M z 
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where we have chosen to keep an energy-independent width. In terms of the physical 
width of the Z into electrons, T e £, and with dLips being the invariant 37 phase space, the 
differential cross section writes as 



+ 



1 d<r(e + e — > 37) 
907r (l67i) 2 dLips 

Pi + 2P 2 _ 
2(3/% + 5/3? 

1 



12tt- 



T% e T 



3" 



z 1 z 

M 2 



Ml 



\Dz(s) 



2(1 _ Xl )(i - X2 )(i - X3 ) - (xlcl (1 _ Xl - ^1) 



1 22 



X 2 



XlX 3 . 



±xlcl (I-X3 



XIX2 . 



1 - xxYii - x 2 y + (1 

-3(l-x 1 )(l-x 2 )(l-x 3 )) 



xi?a 



X3) 2 + (1 



^2) 2 (1 



+ - (x\c{ (1 - 2xt + x 2 2 + x 2 3 ) + x\(? 2 (1 - 2x 2 + x\ + 2$ + x^ (1 - 2x 3 



X3) 



X, 



(3.12) 



Here q is the cosine of the angle between photon i and the beam, q = cos(e 7^) and 
Xi = EJxfs. 



From the above expression we can look at the orientation of the photons relative to 
the beam. For instance, the distribution in the angle ip between the event plane and the 
beam axis, with z = sin(-?/>) (0 < ip < 7r /2) , is given by 



Ida _ 27 g + 2p 2 _ 
tfcU ~ 8 3~K+5pl 



\_#_ 3^ + 14/3i 



9 01 + 2(3 2 _ 



(3.13) 



When specialising to the "fermion-monopole" case, we disagree with ||15|| . In principle it 
is this distribution that can disentangle between different f3. The average value of z may 
be easily computed 

3 15/% + 22/32 



< z >- 



(3.14) 



32 3p 2 + + 5p 2 _ 

The totally integrated cross section at a given energy writes very simply in terms of the 
width measured at the Z resonance: 

37) = ^ T% i ('V \sD z (s)\ 2 . (3.15) 



a NP (e + e- 



M 2 



M 



This shows that, as expected, the cross section grows like s 3 . Hence, at some energy much 
above the Z peak this growth factor can give an enhancement factor of the order of that 
given by the resonance (see below). 
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Of course at the Z peak one has the usual Breit-Wigner formula: 

° (e + e -37) = ^ ^ • (3.16) 

It is also possible to give an exact expression for the s-channel cross section including 
realistic cuts. Besides the cuts on the photon energies and photon-photon separations 
(see Eq. |3.7| ), we can include a cut on the angle, ip, between the event plane and the 
beam to avoid forward-backward events and reduce the purely QED contribution such 
that sin(^) > sin('0o) — z o- The exact analytical formula is lengthy and not very telling. 
With the typical values of the experimental cuts as applied by the LEP collaborations 
the exact formula is well approximated by: 

a -T - 1 - wrw- ( 15(5 * + 9 "- y - 1801 + 2 ^ 6 + ^ o) l 

(3.17) 



4 Monte Carlo Analysis 

4.1 Behaviour of the cross section and effect of the interference 

The above expression can prove handy for a quick estimate of the efficiencies on the 
detection and for the measurement of the Z — > 37 partial width. However, to conduct a 
detailed analysis it is essential to consider the QED 37 background as well as the effect 
of the interference and be able to apply any cut that is needed for a particular detector. 
To include more sophisticated cuts a MC program is necessary. We have written a fast 
MC which includes both the pure QED contribution with the 3 observable photons as well 
as the anomalous terms. The program is based on a calculation of the matrix elements. 
This allows us to study the effect of the interference between the standard part and the 
non-QED part. The very compact expressions are derived through the use of the spinor 
inner product |27|, |3(| |3l| and can be found in the Appendix. By implementing 



the same cuts that we implemented in the analytical formulae we find excellent agreement 
between the analytical results and the output of the MC for all the parts contributing 
to the cross section. We will explicitly show that the interference term can always be 
dropped, hence simplifying the experimental analysis even when one is far away from the 
peak. 

Of course at the Z peak, there is no interference since the tree-level QED amplitude is 
purely real whereas the anomalous is purely imaginary. We find that if the search at the 
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Z peak reveals that the anomalous part is smaller than the QED part (which is the case) 
this implies values for the couplings such that even at LEP2 the interference is negligible. 
Moreover, it is only the term in (3 + that contributes to the interference, the reason being 
that the tree-level cross section for the production of 3 photons with the same helicity is 
known to be exactly zero |2"7|]p|. Since, as remarked in the first section, the classification 
into (3+ and (3- refers to two classes of helicities for the photons, only f3 + interferes. Details 
are found in the Appendix. 

It is always possible to give a parameterisation for the cross section as a function of 
the cms energy, including any cuts. A very simple scaling law parametrisation is obtained 
if the angular cuts are the same at all energies and if the cut on the photon energies are 
implemented through a fixed ratio with the total cms energy. In an obvious notation 
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(4.2) 



For instance, taking the OPAL cuts WA 



cosfl e7 |<0.9 ; # 77 > 20 u 



7 20 



(4.3) 



the values of the cross section crQ,i,+ calculated with Mz = 91.187 GeV, Y z = 2.490 GeV, 
sin 2 9w = Syy = 0.232, are 



(T = 0.705 pb 

o { = 0.082 pb 

a + = 4.181 pb 

a_ = 6.420 pb. 



(4.4) 



Armed with the MC, simulating an experimental environment, we can first exactly 
quantify the effect of the interference. We only consider models with (3 + since there is no 
interference with /3_. As a first illustration we take a value of f3 + /M 4 which corresponds 

"This also holds for same-helicity multi-photon final state processes (29]. 
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Figure 2: Relative effect of the interference for various values of j3 + . Positive f3 + were 
taken. We have preferred, instead of quoting the values of the f3 + , to give the associated 
37 width as calculated from Eq.3.2. 
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to the best limit set on the branching ratio for Z —>■ 37 |13[: Br(Z — > 37) ~ 10 -5 . We 
clearly see (Fig. 0) that around the energies scanned by LEP1, the effect is negligible. 
The relative error introduced by neglecting the interference is below 3% all the way up 
to LEP2 energies. Even for an energy around 300 GeV the effect does not exceed 3%. 
We have also looked at the effect of the interference for a value of (3 + which is 4 times 
smaller (that would give a branching fraction 16 times smaller than the best published 
limit and which could be considered as the ultimate limit from LEP1). The effect up to 
LEP2 energies is below the 1% level. We have also considered values that are 20 times 
smaller; however, foreseeing the accumulated luminosity at LEP1 it would not be possible 
to detect the effect of an anomalous Z — > 37 at LEP. Nonetheless, even for such values 
the effect of the interference is too small even at LEP2 (and beyond) energies. Just to 
make the point we also show that there is negligible effect even when taking larger values 
of the coupling (Figs ||). 

Figures |3] show the merit of sitting at the peak contrasted to the improvement one 
gains by moving to higher (than LEP2) energies, taking for j3 typical values that one may 
hope to measure at LEP1. We see that just above the Z peak the combined cross section 
(QED+anomalous) decreases (due to the Breit-Wigner factor and the 1/s drop of the 
QED cross section) and that it starts bending over and increasing only, unfortunately, 
around the highest LEP2 energies. One conclusion, to which we will come back later, is 
that these figures convey the message that if there is no sign of New Physics affecting 
the Z — > 37 at LEP1 even the highest energies foreseen for LEP2 will not be enough to 
probe this vertex better, even with the higher nominal luminosity of LEP2 (500 pb^ 1 ). 

4.2 Distributions 

The selection of the optimal cuts that should be applied in order to bring out a Z —>■ 37 
effect should be guided by the knowledge of the various distributions in energies and 
angles. These are best given by the MC. We have looked at a few distributions applying 
on both the QED and anomalous parts the same cuts as defined in the previous subsection, 
see Eq. ^L^. We have chosen to show separately the normalised (to one) cross section for 
the anomalous and the QED; the effect of the interference is not taken into account. 
However we have seen that the interference can be safely neglected; moreover, in terms 
of the parameter /?_, this interference is not present, as it is not present at all on the 
Z peak for all operators. Beside the variables we discussed earlier we also looked at the 
distribution in the acollinearity in the transverse plane defined by the most energetic 
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Figure 3: Energy dependence of the cross section for various values of (3 + . The full thick 
curve is the total contribution including the interference while the dotted curve is with no 
interference. Also shown, separately, the irreducible QED cross section. 
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photon (1) and the "medium" photon (2). This is defined from the photon momenta 
transverse to the beam, 2 as 

x = Arcos f SMr-\ . (4.5) 

For all the distributions we show the effect of /3_ and (3 + and one of the other "models" 
defined in the second section in order to see how these models may be disentangled, in 
principle. For the single distribution variables we also show the average value of the 
corresponding variable including the cuts for the QED as well as the anomalous. 

We start with the energy distributions. As argued above there is a clear-cut distinction 
between the QED and the anomalous contribution that is best understood for the least 
energetic photon. In QED e + e~ — > 37 is built up from e + e~ — > 27 where the "added" 
photon can be considered to be a bremstrahl. One should, therefore, expect that among 
the three photons, one is much less energetic than the other two (that will tend to take 
the beam energy). In the case of the Z — > 37, all three photons have on average an 
energy of the same order, thus the energy is shared somehow equally between them. This 
is the reason one of the most salient differences appears in the distribution of the least 
energetic 7. This is well rendered by Fig. [| where two very distinctive spectra stand out. 
This distribution is thus a powerful tool for separating the QED from the anomalous. 
Note, however, that this distribution is almost insensitive to the model of New Physics 
contributing to the s-channel Z —>■ 37 . Therefore based on this remark one could 
devise a "model-independent" optimised cut to unravel the anomalous contributionP]. A 
similar discrimination between the anomalous and the QED occurs in the case of the 
medium-energy photon without "resolving" the models. The distribution in the hardest 
photon shows a less dramatic difference (see Figs. |^). Other distributions in the event 
plane variables relate to the opening angles between any two photons. Although these 
distributions are directly related to the energies (sin 2 (6^/2) = (xj + Xj — l)/xiXj), besides 
bringing out the markedly different structure of the QED and the anomalous, one finds 
that given enough non-QED events a certain discrimination may be perceptible, especially 
in the distribution of the angle between the two softest, Figs. |[ 

Another characteristic of QED events, besides the softness of one of the photons, is the 
collinearity of the photon with the particle that emits that photon. This is not expected 
for the Z —>■ 37 where the photon does not connect to the e + e~ line. Thus another 
**The "knee" that shows up at is purely kinematical and has to do with the ordering of the 



photons. This effect is also manifest in the analytical formula Eq. 3.5 through the occurrence of the step 
function at 23 = 1/2. 
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Figure 4: Normalised distribution in the energy of the softest photon for the case of QED 
as well as for the models defined in the text. Cuts as defined by Eq. |^. j| have been applied. 
The arrows point to the average values of the energy of the softest photon for the models 
considered. 
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Figure 5: As in the previous figure for the medium photon, E 2 and the hardest, E\. 
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Figure 6: As in Fig. [| but in the opening angles between the photons. "3" is the softest 
and "1" the hardest. 
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clear signature is the distribution in the angle between any photon and the beam. The 
QED events, as Figs. [7] show, are clearly peaked in the forward-backward direction even 
after the cuts, whereas the anomalous events are rather central. Note should however 
be made that the softest non-QED photon has also a slight tendency of preferring to be 
in the forward-backward region. In this respect it is really the most energetic photon 
distribution that shows the most marked difference between the QED and the non-QED 
events. Once again, one sees that all models for Z-initiated 37 events give a sensibly 
similar distribution and thus it is highly unlikely that one can resolve any difference 
between the models from the distributions that we have seen up to now. The same can 
be said about the distribution in the acollinearity, once the cuts are applied. 

On the other hand, the angular distribution of the event plane, or equivalently the 
normal to this plane, with respect to the beam axis, does show some interesting structure 
that could be the best way of discriminating between the different models. This is clearly 
seen in Fig. |8| , where we see that not only the average values are well separated but that 
there is a distinct lifting of the "degeneracy" between the models, especially for values of 
the cosine of the angle between the plane and the beam axis above 0.6. 

This lifting of the degeneracy with the devising of more optimal cuts to bring out 
the new physics and attempt to look at its origin can be made much clearer by studying 
distributions in a combination of appropriate variables in the form of scatter plots, for 
instance. In the scatter plot of the energy of the two softest photons, the pure QED 
and the s-channel 37 are confined to two opposite corners (see Fig. |9|). Another obvious 
discrimination is to examine the scatter plot involving the least energetic photon and the 
angle of the most energetic photon with the beam. The scatter plots involving the angle 
of the event plane with the beam show what could, in principle, be the best strategy to 
differentiate between various models. We illustrate this in Fig. [H] by taking as a second 
variable either the energy in the softest or the hardest photon. 

5 Conclusions and Discussion 

We have given a detailed analysis of the most general manifestation of New Physics in 
the coupling of Z — > 37. The analytical formulae should help in quickly estimating the 
acceptances and providing a limit on the operators from the ongoing LEP1 searches. 
We have indicated how to optimise the searches (and the limits) and in the eventuality 
of a signal how to disentangle various effects. Our study also shows that for searches 
around LEP1 energies and with current limits on the anomalous Z — > 37 couplings the 
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Figure 8: Distribution in the cosine of the angle between the beam and the normal to the 
event plane. 




Figure 9: Normalised scatter plots that show the difference between the QED distributions 
and those due to the independent f3+ and {3- couplings. The first is in the energy of the 
softest photon (E 3 ) and the cosine of the angle between the beam and the hardest photon 
(cos9 iz ). The second set of scatter plots is in the energies of the softest (E 3 ) and medium 
photon (Ei). 
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Figure 10: As in the previous scatter plots but as variables: the cosine of the angle between 
the normal to the event plane and the beam, and the energy of the softest or the hardest 
photon. 
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effect of the interference between New Physics and the irreducible QED 37 background 
is abysmally insignificant and can thus be neglected, although it is very straightforward 
to include. 

Let us now ascertain the limits one would derive on the parameters of the effective 
Lagrangian at the end of the LEP1 run. We will assume that each experiment has collected 
a total luminosity J C =200 pb~ 4 and further assume a 100% efficiency on the detection 
of the 37 once the experimental cuts defined by Eq. [D| have been applied. Let us stress 
that since all models, i.e., both f3 + and give sensibly the same distributions and 
that since it will not be possible to reconstruct the helicities of the energetic photons it is 
extremely difficult to distinguish between the different /3's especially if the "anomalous" 
signal is weak. Therefore, the limits will have to be extracted from the total cross section 
after cuts have been applied. It ensues that we will only have access to the combination 
of j3 + and (3- {(3\a + + /^cr_) where cr± are defined in Eq. |4J] and depend (mildly) on 
the cuts applied. Without the cuts one has the combination that appears in the width, 
(3/3++5/3_). With a healthy statistics and a good signal it could be possible by judiciously 
varying the cuts to "reconstruct" the correct combination of (3- and j3 + . In what follows 
we will restrict the discussion to one parameter. Requiring a 5a deviation, with only the 
statistical error taken into account, one finds for f3_ for instance 

4a 2 



m 4 



< 0.215 (0.256 for y £ = 100 pb" 1 ). (5.1) 



The increase in luminosity does not, unfortunately, tremendously improve the limits, as 
the effect of the New Physics is quadratic in the couplings and therefore the sensitivity 
scales only as £ 1//4 . This remains true, as we have seen, even when one moves away from 
the Z peak. Thus, one may wonder, since the QED cross section decreases with energy 
as 1/s while the anomalous grows as s 3 , whether one could set a better limit at LEP2. 
We have already answered by the negative. We show in Fig. |ll| how the significances 
are changed as the energy is increased while keeping the same cuts and normalising to a 
common integrated luminosity of £ =100 pb -1 ]^. As the figures demonstrate, for various 
values of (3 + (in the range set by LEP1), the statistical significances are never better at 
LEP2 than what they are at LEP1. One will take full advantage of the energy increase 
of the anomalous cross section, bettering the resonance enhancement, only for energies 
around 300 GeV. Taking again the operator defined by /?_, and with LEP2 meaning 

TTOf course, as the energy increases one expects the luminosity to improve as well, however we stress 
again that since the effect of the interference is marginal, the significance can be easily calculated through 
the scaling law £ 1//4 . 
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■s/s = 190 GeV, C =500 pb 1 a limit based on a 5a (statistical only) deviation gives 

4a 2 ~ 

< 0.365, (5.2) 



m 4 



ie., worse than at LEP1. However, at a Next Linear Collider (NLC) operating at 350 GeV 
with £ =10 fb _1 we find that the LEP1 limit is improved by as much as an order of 
magnitude: 

4a 2 ~ 

—rW\< 0.025. (5.3) 

m 4 

This said, it is not excluded that the strength of the 4-7 vertex is much larger than 
that of the Z — > 37 vertex, in which case 37 production, proceeding through photon 
exchange, will be much better studied than at LEP1, since it does not receive any res- 
onant enhancement. Therefore we should urge the LEP2 collaborations to scrutinise 37 
production. The topology of these 37 events is exactly the same as those originating 
from the e + e~ — > Z —>■ 37 that we have studied here; in particular, all the normalised 
distributions are the same. To adapt what we have investigated in this paper to the effect 
of a 4-7 vertex in 3-7 production one simply makes the replacement j3± — > (3± together 
with the change of the Zee vertex into the electromagnetic one and the replacement of 
the Z propagator. For instance assuming that (3± = sw/cwP± as given by a model with 
a strongly interacting SU(2) singlet fermion, then the 4-7 strength is larger than that of 
the Z — > 37 . For this particular case the combined effect of the photon and the Z is 
accounted for by the factor 



V 9l s J \ s w s 

9l + 9r 



(5.4) 



that multiplies the results of taking into account the Z-channel only. At 190 GeV this 
enhances the Z cross section by almost seven fold. At LEP2, there is also the possibility 
of checking for the 4-7 vertex through 27 production in the (peripheral) 77 processes. A 



dedicated investigation of this coupling for LEP2 and beyond is under investigation ||32 |. 
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Figure 11: Statistical significances as a function of the cms energy assuming a common 
effective luminosity L =100 pb^ 1 . Three typical values of the branching ratio are assumed 
with f3+ > 0. 
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Appendix: 



A Helicity amplitudes with the spinor inner product 
technique. 

With the cuts that we had to implement, the helicity amplitudes could be evaluated by 
neglecting the electron mass. There is a variety of methods for the evaluation of the 
helicity amplitudes |27], PEfl . The most powerful are those that have been developed for 
the calculation of multiparticle processes especially in massless QCD [[!(], [H], |28| and are 
thus well adapted for the process at hand. The amplitudes are written in terms of Spinor 
Inner Products. We have performed all the helicity amplitudes calculation with the help 
of FORM. We have also checked our results against known results for the QED part 
and against the usual trace technique summation for the matrix element squared. For 
the helicity technique we have made full use of gauge invariance by choosing the most 
appropriate choice of gauge for the polarisation vectors to render the expressions as simple 
as possible. This is a huge gain on the usual squaring technique. The helicity amplitudes 
can all be expressed though the inner product S(p, q), where p, q are (light-like) momenta. 
Following the handy notation of |31|, the massless spinors are written as 



\p± >= u (p, ±-) = v(p, =f- 
<p±\=u(p,±-)=v(p,^-) (A.l) 



then 



S(p,q) = < p - \q+ >=< q + \p- >*= -S(q,p) ; 
\S(q,p)\ 2 = 2p-q ■ (p 2 = g 2 = 0). (A.2) 

S is antisymmetric, unlike in [RDl . 

We decompose the amplitudes into the standard QED part and the two non-interfering 

j3± as 



The helicity components are labelled as .4.oi,<72;Ai,A2,a 3 > where the first labels refer 
to the helicity of the e~ and e + . Chirality (m e = 0) only allows the helicity ampli- 
tudes A a - a -\ lt \ 2t \ 3 . One then only needs the expressions for A+-- t ± t ± t ±, A + -- i+ -- and 
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A + .-■++-■ For QED we recover the known result that the amplitude with the like-sign 
photon helicities vanishes: 

•4?SU ± = (A.4) 
while the two others are simply given by 

i qed _ S(p 1 ,p 2 )S(p 1 ,k 3 ) 2 



A) 



S(pi, k^Sfa, k 2 )S(p 2 , h)S(p2, k 2 ) 



aQed ( Sjpup^Sjp^ks) 2 \ 

" \s{p 1 M)sipiM)sip2M)sip2M)) 1 ' 

For the anomalous part, B~ only contributes to the like-sign photon helicity amplitude 
whereas B + only for mixed helicity states: 

Bt- i+++ = K-,— = = = o (a.6) 

and 

BX---+ = 9r s( Pl ,p 2 y(s(k 3 , Pl )s(k u k 2 y) 2 

= -g R SfaptXSfanySfakd) 2 (A.7) 

$-+-,+++ = g R S( Pl , P2 y (S(k 2 ,k 3 ) 2 S(k l7Pl ) 2 

+S(h, h) 2 S(k 2 , Pl ) 2 + S(h, k 2 ) 2 S(h, Pl ) 2 ) 

= -9r s(pup 2 ) (s(k 2: k 3 ) 2 s(k uP2 ) 2 

+S(h, k 3 ) 2 S(k 2 ,p 2 ) 2 + S(h, k 2 ) 2 S(k 3 ,p 2 ) 2 )* . 

(A.8) 

All other amplitudes are found by permutation in the photons and parity conjugation 
that reverses all signs of the helicity. For the Z mediated amplitudes this conjugation 
should also be accompanied by allowing for g R — > g L . 

The numerical generation of S(p, q) is done through the definition 

Sip, q) = - 7 1 + yj (go + qaj - -7= + = \J {Po + Ps) ; P = (po,Pi,P2,Pa)- ( A -9) 
(po + P3) J (go + 93) 
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Where to find the generator? 

All the cross sections and distributions found in this paper have been obtained by numer- 
ical integration of helicity-amplitude-based matrix elements. This was necessary in order 
to show that the interference between the QED background and the anomalous signal is 
negligible for energies up to LEP2 and for reasonable values of the couplings. However, 
from an experimentalist point of view, this is not exactly what is needed. To conduct a 
detailed experimental analysis, it is more important to have an event generator for the 
QED background and a second one for the signal. The former can already be found on the 
market [j33| and we now provide a new generator for the anomalous signal. It can be found 
at Iapphp0.in2p3.fr/pub/preprints-theorie/ee3gammagenerator.uu using anonymous ftp. 
The uu-encoded file contains five files of which only one is important. The four others are 
there only to run a small demonstration program to see if it runs well on your system. 
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